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^^ \ Nuclear modifications to the Drell-Yan dilepton production cross sections in p+A 

1-^ I and A + A collisions in the leading twist approximation are caused by nuclear 

effects in the parton distributions of bound nucleons. For non-isoscalar nuclei, 

Qh. isospin corrections must also be considered. We calculate these effects for \i+A 

(— I \ and Pb+Pb collisions at CERN SPS energies. Our goal is to place constraints 

on nuclear effects in sea quark distributions in the region x^Q.2. We show that 

the net nuclear corrections remain small for 'p+A collisions at ii^iab = 450 GeV. 

rS \ However, in Pb+Pb collisions at E'lab = 158y4GeV, effects of ^ 20% are predicted 

C^ ' at large M. The data collected by the NA50 collaboration could thus be used 

to constrain the nuclear effects in the sea quark distributions in the region of the 

EMC effect, x>0.3. 
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1 Introduction 

In this work, we study perturbatively calculable QCD processes in high energy p+A, 
A + A and A + B collisions. In the search for the Quark-Gluon Plasma (QGP) in 
ultrarelativistic heavy ion collisions, these inclusive hard processes provide reference 
cross sections for detecting specific signatures of the QGP, such as suppression of heavy 
quarkonia, production of thermal dileptons and photons, and energy losses of jets. 

At high energies the inclusive differential cross sections of hard processes in nuclear 
collisions can be computed (in leading twist approximation) from a factorized form 
consisting of nuclear parton distributions and partonic cross sections. At this level 
all nuclear effects are contained in the nuclear parton distributions, which obey the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations of perturba- 
tive QCD (pQCD) ^. At higher orders in pQCD, absorption of (1/e) singularities into 
definitions of parton distributions is scheme- dependent (usually MS) , which retains the 
same universality properties as in the case of hard processes for free nucleon scatterings. 
Consequently, the same nuclear parton distributions can be used to compute different 
hard scattering cross sections in nuclear collisions. This is the obvious motivation to per- 
form a consistent DGLAP analysis of nuclear parton distributions, as done in 0, |^, ^. 

Symbolically, the inclusive hard scattering cross sections for producing a particle c in 
a collision of nuclei A and B can be written as 



da{Q^, ^/s)AB^c+X = Yl 



i,3=q,q,9 



ZAff'\x,,Q') + {A-ZA)f:'\x,,Q') 



ZBff''{x2, g2) + {B- ZB)ff''{x2, Q2)] ® da{Q^, Xi, X2),,-.c+x (1) 

where cr^Q"^, Xi, X2)ij^c+x is the perturbatively calculable cross section at a large momen- 
tum (or mass) scale Q ^ ^qcd ~ 200 MeV, Xi^2 ~ Q/V^ are the fractional momenta, 
/f (xi, Q^) is the distribution of parton species i in a proton (neutron) of the nucleus 
A, and correspondingly /J is that for the nucleus B. The number of protons in A{B) 
is denoted by Za{Zb)- For isoscalar nuclei, the parton distributions of bound neutrons 
are obtained through isospin symmetry (as in the case of unbound nucleons), f^L. = f^U 

and /^As = f^Ly This is expected to be a good approximation for non-isoscalar nuclei 
as well. Therefore, one may formulate the studies of nuclear parton distributions simply 
in terms of those in bound protons, which are denoted here as /f = fi/A- It is con- 
venient to define the nuclear effects in parton distributions in terms of the ratio of the 
distribution of the parton species i in a bound proton to that in a free proton, 

'^ '^^- Mx,Q')- ^^ 



Information on the nuclear parton distributions is primarily obtained by deep inelastic 
lepton-nucleus scattering (DIS) 0-|]l3l and by Drell-Yan (DY) dilepton production in 



proton-nucleus collisions |jT^, |T^. The nuclear structure functions F2{x,Q'^)^ measured 
in DIS P]-|jl3[ differ from those in free nucleons. The ratios R^^ = ^Fg^^/^F^, where 



deuterium D approximates an average free nucleon, show clear and systematic devia- 
tions from unity in various regions of Bjorken-x: shadowing (i?^^ < 1) at a; ^0.1, anti- 
shadowing (i?^ > 1) at 0.1<x<0.3, EMC effect (/?4 < 1) at 0.3 < a; < 0.7, and Fermi 
motion (i?^^ > 1) as s ^ 1 and beyond. The New Muon Collaboration (NMC) high- 
precision measurements of the F2 structure function ratios for tin vs. carbon, F2'^ / F2 
have also revealed a dependence on the virtuality scale Q^ at small values of x. These 
measured modifications of nuclear structure functions directly imply modifications of 
parton distributions in bound nucleons. 

Just as in the QCD-improved parton model (in lowest order) the structure functions 
can be written in terms of parton distributions, 

Fi{x. Q') = E el{Z[xf^/^ix, Q^)+xf?^^{x, Q')] + iA-Z)[xf-/^ix, Q^)+xf^'^{x, Q^)]}, 

(3) 
where q is the quark flavour and e^ is the corresponding charge. In the DGLAP analysis 

of nuclear parton distributions 0, |^, || it is assumed that the distributions /f are 
factorizable at a sufficiently large initial scale, Qq ^ ^qcd- Once the input distributions 
are given at Qq and at x > Xmin-, their evolution is predicted by the DGLAP equations 
aX Q > Qq and x > Xmin- As a result, the ratios Rf{x, Q"^) depend both on x and Q^. In 
analogy with the global analyses of the free parton distributions, the key problem is then 
to determine the (nonperturbative) initial distributions /f (x,(5o)- To constrain these, 
further information is needed. This is provided by the DIS measurements mentioned 
above (NMC ||-||, SLAC ^, E665 fll], g) and by the Drell-Yan data from the 



E772 and E866 collaborations in p+A collisions |n, 15|. In addition, conservation of 
momentum and baryon number serve as further constraints. We emphasize that the 
measured Q^ dependence of the ratio F2^ / F2 P is also reproduced very well by the 
DGLAP evolution §. 

In the DGLAP analysis 0] of the nuclear parton distributions the nuclear effects 
were expressed in terms of free parton distributions which were assumed to be known, 
i.e. obtained from a set of distributions such as CTEQ, GRV, MRS etc. The absolute 
distributions from different sets of free parton densities may differ by a fairly large factor^ 
and, consequently, these differencies will be reflected in the absolute nuclear parton 
distributions as well. The ratios -R^(x, Q^), however, vary only within a few percent 
from set to set, as shown in [Q. Therefore, for computing hard processes in nuclear 
collisions with nuclear effects in the parton distributions, it is a good approximation to 



^ideally of course, there would be only one best set {fi{x, Q^)} 



use universal ratios Rf{x, Q^) which are independent of specific free parton densities. A 
parametrization of Rf{x,Q^), "EKS98", was prepared in Q] for general use, and it is 



available in |16] and now also in the the CERN PDFLIB library of parton densities [17 



Some uncertainties, however, remain in the determination of the nonperturbative 
input distributions at Q^. In this paper, we will focus on constraining uncertainties in 
the input sea quark distributions in the region a; ^0.2. .0.3, i.e. approaching the region 
of the EMC effect, where the ratios Rp^ measured in DIS are dominated by valence 
quarks. Our goal here is to study to what extent the NA50 Drell-Yan data for p+p and 
p+A collisions at E'lab = 450 GeV and Pb+Pb collisions at i^iab = 158y4GeV at the 



CERN-SPS |]T8[ can be used to constrain the EMC effect for the input distributions of 



the nuclear sea. 

2 General properties of nuclear corrections 

Let us first discuss in some detail how available data and sum rules constrain the in- 
put nuclear parton distributions, or equivalently the nuclear modifications Rf{x, Ql), in 
different regions of x in the DGLAP approach [^. Since the DGLAP analysis a pertur- 
bative, the scale evolution must be limited to the region Q > Qo ~ 1 GeV.[] Some hints 
are given, however, by DIS measurements in the non-perturbative region Q < Qq. As 
illustrated in Fig. 1 of [Q, the experimental constraints from DIS and DY are not given 
along a fixed value of Q^ - as would be preferable for the DGLAP initial conditions - 
but in certain kinematically correlated regions of x and Q^. Furthermore the data from 
DIS and DY are typically in distinct kinematical regions. For these reasons the input 
distributions must be constrained by using a recursive procedure similar to the global 



analyses of free parton distributions [23, 25 



In first approximation the input nuclear effects for valence and sea quarks can be 
assumed to be separately flavor- independent: -R^^(a;,<5o) ~ ^dy^-^^Qo) ~ Rvi^^Qo)^ 
and R^{x,Ql) ^ Rf{x,Ql) ^ Rf{x,Ql) = Rs{x,Ql) 0. Thus only three independent 
input ratios must be constrained at Qq: Ry, Rg and Rq. 

Quarks and antiquarks 

• At a; <; 0.3 the structure function F^ is dominated by valence quark distributions. 
The DIS data for Rp^ therefore only constrains the magnitude of the EMC effect 
and the Fermi-motion in Ry but not those in Rg or in Rq. 

• At 0.04 < X < 0.3 the DIS and DY data give bounds for R^ and Ry but in different 
regions of Q^, (see Fig. 1 of [H]). 



^In Q Qo = 1-5 GeV was chosen. 



• At a;<0.04 there are DIS data for the ratio Rp^ available down to a; ~ 5 ■ 10~^ 
in the region Q^l GeV relevant for the DGLAP analysis. Conservation of baryon 
number forces the nuclear valence quarks to be less shadowed than the sea quarks. 

• In the DIS data for Rp^ at x ^ 5 ■ 10^^ one enters the nonperturbative region Q ~ 1 
GeV. A saturation behaviour of Rp^ in x — ;> is observed along the experimentally 
probed values of Q^ ^, |Tl| . Provided that the sign of the slope of the Q^ depen- 
dence of Rp^ in the nonperturbative region remains the same (positive) as what 
is measured at x ~ 0.01 in the perturbative region 0, a saturation behaviour, 
i.e. a weak dependence of Rp^ on x, can also be expected at Q^. Constraints are 
then given by the DIS data in the non-perturbative region, in the sense that the 
data give a lower bound for Rp^{x,Qq). Since at small values of x the sea quark 
distributions dominate over the valence distributions, Rg is also constrained by the 
DIS data while the shadowing in Ry is restricted by baryon number conservation. 

Gluons and sea quarks 

• The scale dependence of Rp^ at small values of x is directly connected with shad- 
owing of gluons: the more deeply gluons are shadowed, the slower is the evolution 
of Rp^. The ratio Rq, can thus be constrained by the measured Q^ dependence of 
F^ as done in |1^, |^. Since the Q^ dependence is not very strong, high-precision 
data is needed. In practice only the NMC data for Q^ evolution of F2^ / F2 can 
be used for constraining the input nuclear gluons at 0.02 <x< 0.2. 

• At x<0.02 it can be assumed that Rq^x.QD ^ Rp^{x,Ql) for x -C 1. This 
remains true within about 5 % even after the DGLAP evolution from Qo ~ 1 GeV 
to Q~ 100 GeV i. 

• At x^O.2 there are no direct experimental constraints available but conservation 
of momentum together with arguments for stable evolution can be used. In this 
region one is sensitive to the small tails of the gluon distributions, and the existence 
of the EMC effect for nuclear gluons cannot be deduced based on the momentum 
sum rule alone. The evolution of gluon distributions are, however, affected by the 
valence quark distributions (but not vice versa), so since an EMC effect exists 
for the valence quarks, one will be generated for the gluons as well. The evolution 
equations for gluons and sea quarks are mutually coupled, so an EMC-like depletion 
will be generated for the sea quarks through the DGLAP evolution p[. Then if 
the nuclear ratios Rf do not move away from their input values very rapidly, it 
is a plausible first approximation to have an EMC effect both for the input gluon 
and sea quark distributions. For the sea quark modifications, which is the major 
subject of this paper, a simple assumption of Rs{x^ 0.3, Ql) ~ /2y(x<; 0.3, Qg) 



was made in for the input modifications of the sea quarks. We now move on to 
study the effects of this approximation on the Drell-Yan dilepton cross sections in 
nuclear colhsions at SPS energies. 

3 Nuclear effects in DY production 

Nuclear effects in Drell-Yan dilepton production can be divided into two classes: first, 
there are "genuine" nuclear effects arising from the dynamics of the nuclear collision. 
These include nuclear modifications of parton distributions. Second, even without any 
nuclear effects in the parton densities or in the collision dynamics, the DY cross sections 
in p+A collisions (normalized per A) differ between isoscalar and non-isoscalar nuclei 
due to the different relative numbers of protons and neutrons. It is often convenient to 
use deuterium D as a reference, since it is approximately a sum of a free proton and a free 
neutron. Then any observed deviation of the DY cross section for p+A in isoscalar nuclei 



1^ 1 from that for p+D can be interpreted directly as a genuine nuclear effect, such as a 
nuclear modification in the parton distributions. For non-isoscalar nuclei, ratios of DY 
cross sections for p+A to those for p+D or p+p always show additional isospin effects. 
For the purposes of comparison of the DY cross sections from p+A with those from 
p-|-D and p-(-p, we write the isospin symmetric part of the parton distributions of the 
nucleons in a nucleus with Z protons separately, 

Z/f/^ + {A- Z)ff^ = ^(/f/^ + /f/^) + {Z- |)(/f/^ - /f/^). (4) 

The inclusive cross section for the production of the Drell-Yan dilepton pairs of invariant 
mass M and rapidity y in p+A collisions can then be written in the lowest order as 

q=u,d 

where g^ = fq{xi, Q"^) with the momentu 

Q^ = M^. The ratio of the inclusive Drell-Yan cross section in a p+A collision vs. that 

in a p+p or p+D collision now becomes 

"""^^ ^ ' ^dafy/dMHy 

= {4[mi(m^ + d^) + ui{u^ + 4)] + [di{d^ + u^) + (Ji(4 + u^)] + 4si4 + ■■■}/NEy 

97 

+ (— - l){A[u,{4 - d^) + u,{u^ - 4)] + K(J^ - u^) + Ji(4 - u^)]}/NEy (6) 



2 


'qM"' + f2'^)+<iMl"' + q2'^) 


+ 


'qM'^~q2'^)+qi{ql'^-q2"') 


}, (5) 


m 


fractions Xi^2 = -%e=''^ and a 


scale choice 



where the denominator N^y only contains free parton densities. For deuterium [B = 
D = 2) NEy is 

N£y = 4[mi(m2 + J2) + mi(m2 + d2)] + [di{d2 + U2) + di{d2 + M2)] + 4siS2 + ... (7) 
and, correspondingly, for the proton {B = p = 1) it is 

Ndy = NEy + H 4[mi(m2 - d2) + mi(m2 - d2)] + [di{d2 - M2) + diid2 - U2)]. (8) 

q=u,d 

It is evident that for isoscalar nuclei, A = 2Z, the ratio Rdy is unity in the absence of 
nuclear modifications in the parton densities. At large rapidities (large xp) X2 <C xi and 
the ratio RjJy is sensitive mainly to the nuclear effects in the sea quark distributions. In 
the following, however, we are interested in the central rapidities y ~ 0, so Xi ~ X2 and 
the ratio Rjyy thus reflects the nuclear effects both in the sea quark and in the valence 
quark distributions. 

For non-isoscalar nuclei on the other hand, even without nuclear modifications in the 
parton distribution functions, the isospin corrections ~ (^ — 1) must be considered. In 
fact, they depend quite strongly on the specific set of parton distributions used in the 
calculation. In the oldest sets, u = d was assumed. Since A > 2Z, this would lead to 
RjJy (x, Q^) < 1 in the absence of nuclear effects in the parton distributions. However, 
the ratio u/d is experimentally different from unity: NA51 Collaboration measured this 



ratio for the first time, reporting a value of u/d = 0.51 ± 0.04 ± 0.05 at x = 0.18 pO 
This fact was taken into account in the subsequent sets of parton distributions, such as 
MRS94 0, GRV94 [|2[ and CTEQ94 pg. Later on, E866 Collaboration at Fermilab 



measured this ratio with higher accuracy [24|. The modern parton distribution sets now 



include the u/d asymmetry. As a result, the behaviour of the isospin corrections for 
the Drell-Yan process calculated by using the recent parton distribution sets differ from 
those calculated by using the old ones. To demonstrate this, we plot in Fig. |I| the ratio 
Rdy computed with the parton densities of the free proton. For tungsten, A = 184 and 
Z = 74. The ratio is shown as a function of mass M at i/i = 30 GeV and ycm = 0.04 
for four different sets of parton distributions. 

Next, we study how the nuclear effects in the parton distributions and the isospin 

effects are reflected in the ratios -^ d^M /\ 'dM ^^'^ \ Im I 'mi • These ratios can be 
formed from NA50 experimental results. They have measured inclusive dilepton produc- 
tion in p+p, p+D, p+ |Be and p+^^|W colhsions at Eiab = 450 GeV (a/s = 30 GeV) at 
the CERN SPS in the rapidity range 3 < yiab < 4 (-0.46 < y^m < 0.54) and mass M 
around the J/\E' peak. For M > 4 GeV, the mass spectrum is dominated by Drell-Yan 
dileptons. 

For the following calculations, we integrate cross sections (|^) over the NA50 rapidity 
bin, and form the above ratios. For the free parton distributions we use the set MRST98 
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A=184, Z=74, v^s=30 GeV, ycm=0-04 
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Figure 1: Isospin effects in the ratio of DY cross sections -R^iy of Eq. (^) plotted as a function 
of mass M for ^f^W at ^/s = 30 GeV and ycm = 0.04, computed with parton distributions from 
MRST98 (central gluon) l2| (sohd), CTEQ4L Q (dotted-dashed), GRV94LO ^ (dotted) 
and GRVLO [gg] (dashed). Nuclear effects in the parton distributions are not included. 



(central gluons) [^ and the nuclear effects in the parton distributions are taken from 
the EKS98 parametrization [Q. 

In Fig. these ratios are plotted in lowest order for ^Be (upper panel; solid, dashed) 
and ^ylW (lower panel; solid, dashed). The analysis for the nuclear effects i?f(x, Q^) 
is only a leading order one, so strictly speaking it should be used only together with 
the leading order parton densities. The Q^ evolution of the ratios Rf{x,Q'^), however, 
is relatively slow, so the ratios given by EKS98 also serve as a first approximation for 
the nuclear effects in the next-to-leading-order (NLO) parton distributions. Keeping 
this source of uncertainty in mind, we have also computed the Drell-Yan cross sections 
in NLC)0 |28|. The ratios resulting from the NLO computation are also shown in Fig. 
H for 4Be (upper panel; dotted-dashed, dotted) and ^74W (lower panel; dotted-dashed, 
dotted). As seen in the figure, the LO ratios are a good approximation to the NLO 
ratios. 

In NLO, 0{a'^as), the Drell-Yan cross section consists of quark-antiquark annihila- 
tions with an emission of one real gluon; qqA,qqA —>■ 91*, and gluon-initiated Comp- 
ton scatterings; qgA,gqA -^ Ql* and qgA,gqA -^ Ql*, and one-loop corrected quark- 
antiquark annihilations interfered with the LO annihilation. The total NLO contribution 



^Fortran code from P.J. Rijken and W.L. van Neerven is used. 
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Figure 2: The M dependence of the ratios of Drell-Yan cross sections \-jjj- for p+A-collisions 
at y/s = 30 GeV and 3 < yiab < 4. Upper paneh p+Be over p+p for LO (sohd) and for NLO 
(dotted-dashed), p+Be over p+D for LO (dashed) and for NLO (dotted). Lower panel is the 
same for p+W. Shadowing and isospin corrections are taken into account 



to da/dM'^dy can be written as a sum da^io = <^s+v + o'Hgq + o'gq, where <7s+v are the 
virtual corrections summed together with the soft emission of a gluon in the qq annihi- 
lations, aHqq contains the emissions of hard gluons in qq annihilations, and agq accounts 



for the Compton processes |]28|. In the kinematical range studied here, (Js+v dominates 
the net NLO contribution. The sign of anqq changes: anqq/crs+v ~ +0.04... — 0.39 
for M = 10... 1.5 GeV and ?/cm = 0.04. The Compton term agg remains negative and 
(Tgq/crs+v ~ —0.12... — 0.09. Thus the net NLO effect is always a sum of partly canceling 
contributions and Kdy = ^Md / dlid ~ 1-6. ..1.7. In the ratios Rjjy and Rj^/y, however. 



1.04 



23456789 10 
M [GeV] 

Figure 3: The ratio of shadowed over non-shadowed Drell-Yan cross sections ^j^ for p+Be 
LO (solid), p+Be NLO (dotted-dashed), p+W LO (dashed) and p+W NLO (dotted) colhsions 
at ^ = 30 GeV and 3 < yiab < 4. 

the i^-factors largely cancel and the ratios remain very close to those computed in the 
LO, as seen in Fig. 0. 

The ratios in Fig. contain both the nuclear effects in the parton densities and 
the isospin corrections. To see the effect of the nuclear parton distributions alone, we 

have plotted in Fig. ^ the ratio of the cross sections \-^^ computed with and without 
nuclear effects in the parton distributions. Both the LO and NLO ratios are shown. 
The isospin corrections are taken into account in all cross sections; thus without nuclear 
modifications of the parton densities all the ratios in the figure would reduce to unity. 
By comparison with Fig. ^, and from Fig. |^, we conclude that the isospin corrections to 
the ratios Rjjy remain small in magnitude when the MRS98 distributions are used. 

Figs. |] and |] indicate that the net nuclear effects in DY caused by the nuclear 
modifications of parton densities are not very dramatic in p+A collisions at the SPS 
energy i?iab = 450 GeV/c in the kinematic region 1 GeV ^ M < 10 GeV, 3 < |/iab < 4. 
There are two reasons for this: first, at the corresponding values of X2 the nuclear effects 
for the sea quarks are small. Second, since one is predominantly in the antishadowing 
region for valence quarks, there is a cancellation in the net nuclear effects in the ratio 
shown. To demonstrate this, and to show to what extent the ratio pA/pD reflects the 
nuclear modifications of the sea and valence quarks, we show in Fig. H the nuclear effects 
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Figure 4: The ratios R^{x2,M'^) (solid) and R^^{x2,M'^) (dotted) for if|W as a function 
of invariant mass M (lower x-axis) and X2 = — :=e~^''™ (upper x-axis) for ^/s = 30 GeV and 
ycm = 0.04. The ratio of shadowed vs. non-shadowed LO DY cross sections (dashed) is the 
same as the dashed curve in Fig. y. 



R^^{x2, M^) and R^{x2, M^) as a function of mass M, with X2 computed in the middle 

"^ where ?/cm = 0.04. The net effect 






of the accepted rapidity bin (LO only), X2 

in the DY ratio pW/pD is shown again by the dashed curve. The values of X2 can be 
read off from the top of the figure. Note here that the scale evolution of the sea quark 
modifications is taken into account. 

To constrain the nuclear effects of the sea quarks at larger values of x, we must consider 
lower energies. The NA50 Collaboration at CERN has measured dilepton production 
in Pb+Pb collisions at Ei^b = 158 AGeV {^/s = 17.2 GeV). This offers us a better 
chance for constraining the EMC effect in the input sea quark distributions. To illustrate 
the sensitivity of the DY dilepton cross sections da^^'^^^ / dM"^ (LO, integrated over 
3 < y\ah < 4) to the assumption of the EMC effect in the input modifications of the 
nuclear sea quarks, the ratios 



dcrT'" 



daP^+P^/dM^ 



(9) 



are plotted in Fig. ^. The cross section 



dtjPb + Pb 



is the DY dilepton cross section which 



includes the isospin effects as in Eq. |l] but no nuclear effects of the parton distributions. 



11 



Pb+Pb 



The cross section — J^j2 — similarly includes the isospin effects but is computed by using 
three different scenarios for the nuclear effects of the sea quarks. First we take into 
account all the nuclear effects as given by the EKS98. The resulting ratio is shown by 
the solid curve in Fig. |^. Second, we keep the nuclear effects of the valence quarks in ac- 
cordance with EKS98 but relax the assumption of the EMC-effect in the sea component. 
We take the ratios Rf{= Rf^) from EKS98 at x < 0.1 (where these ratios are constrained 
by experimental data) but interpolate R^ linearly from x = 0.3 to the region of Fermi 
motion x > 0.8 without an EMC effect]^. The resulting ratio is shown by the dashed curve 
in Fig. ^. Finally, we switch off all nuclear effects in the sea quarks and antiquarks by 
setting R^ = 1. We emphasize that the last scenario is actually unphysical as it violates 
the available constraints, and that it is meant only for comparison purposes, to see the 
effects of the nuclear corrections in valence quarks alone. The corresponding ratio is 
shown by the dotted curve. Note that the deviation from unity of the ratios in Fig. |^ 
directly shows the effects of the nuclear modifications in the parton distributions. The 
net effect of nuclear parton distributions is now clearly larger than in the p+A case. 

In principle it should be possible to form this ratio from the measurements, by taking 
the numerator directly from Pb+Pb data and the denominator from p+p and p+D 
data. To our knowledge, however, no experimental data presently exists for Drell-Yan 
production in p+p or p+D collisions at _Eiab = 158 AGeV/c. Thus, some additional input 
is needed in order to form the experimental ratio (§). One possibility is to compare the 
measured DY cross section in Pb+Pb directly with a purely theoretical calculation. 
Alternatively, the denominator in Eq. (^ could be formed from other p+p or p+D 
data (e.g. data from NA51 collaboration at £^iab = 450 GeV/c ||20|) corrected to i?iab = 
158 AGeV/c, based on theoretical cross sections without nuclear effects in the parton 
distributions. In both cases modern parton distributions, in which constraints from the 
measured ratio u/d are included, must to be used in order to correctly include isospin 
corrections. The main uncertainty would be the overall normalization, which is common 
to both cases. As shown by Fig. |^, the experimental ratio is expected to be unity within 
a ~ 5 % uncertainty at M ~ 4 GeV. If a direct measurement could be made at a single 
M-value, it could be used to fix the overall normalization of the ratio. Then the slope 
of the ratio towards larger values of M would give direct information of the EMC effect 
in the sea quark distributions. 

Our study shows that the expected nuclear effects are ^ 20 % at masses larger than 7 
GeV. Relevant constraints would be obtained for the EMC effect of the sea quarks and 
antiquarks if the precision is ~ 10 %. 



^This procedure would cause some inconsistency with the EKS98-niodification of gluons due to the 
scale evolution of the parton densities but as we now do the DY computation in the LO only, gluons 
are not involved. 
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Figure 5: The ratio @ of the distributions da/dM"^ at 3 < yiab < 4 computed with and 
without the nuclear effects in parton distributions for Pb+Pb collisions at E'lab = 158 AGeV. 
Solid curve: Rq taken from the EKS98 parametrization. Dashed curve: Rq interpolated over 
the region 0.1 < x < 0.85 by assuming no EMC effect. Dotted curve: i?^ = 1. Isospin effects 
are included. The values of X2 = — ^e"^''™ and xi = — ^e^''™, computed in the middle of the y 
bin can be read off from the upper x axis. 

4 Conclusions 

We have studied the sensitivity of the Drell-Yan cross sections in SPS nuclear collisions to 

isospin corrections and to the nuclear modifications of the parton densities in the regions 

of antishadowing and EMC effect. Our aim has been to find a way to constrain the nuclear 

modifications of the sea quark distributions in order to improve the determination of the 

input modifications in the DGLAP analysis [Q of nuclear parton distributions. We have 

shown that the Drell-Yan dilepton data collected by the NA50 collaboration at CERN 

SPS in Pb+Pb collisions at ii^iab = 158 AGeV/c would be suitable for constraining 

the EMC effect in the input distributions of the sea quarks provided that a sufficient 

^^pb+Pb Pb+Pb 
precision, ~ 5. ..10% is reached in forming the ratio J^j2 / dM'^ at M>5 GeV. 

We have also shown that for p+A collisions at i^iab = 450 GeV, in the kinematic range 

3 < l/iab < 4, the net effects due to nuclear modifications of the parton densities are small. 

This is because typical values of x remain in the region where the nuclear effects in the sea 

quarks and valence quarks largely cancel. Even at the highest masses studied, M ~ 10, 

where the typical x for E'lab = 450 GeV is ~ 0.3, the nuclear parton distributions modify 
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the p+W cross sections only by 6 %. This sets the minimum precision required for such 
an experiment to constrain the nuclear sea quark distributions in the antishadowing 
region. We have also shown that the isospin effects are small, provided that modern 
parton distributions, where u ^ d, are used. 

Based on the data from p+A collisions, it is often assumed that Drell-Yan behaves as 
A^ in nuclear collisions. From our results for the SPS, Figs. ^ and ^ we see that strictly 
speaking this is not the case but the deviations remain fairly small, and within a 5 % 
uncertainty in the cross sections the deviations can be neglected. However, at higher 
energies, such as i^iab = 800 GeV in the Fermilab E772 experiments, the shadowing cor- 
rections at smaller values of x become important and have been experimentally observed 

3,0. 



Finally, let us comment on the consequences of these nuclear effects in the analysis 
of the Drell-Yan cross sections in Pb+Pb collisions measured by the NA50 collaboration 
in connection with J/\E' suppression [Q. We have shown that the slope of the invari- 
ant mass distributions of the DY pairs is affected by nuclear effects: Fig. ^ indicates 
that corrections of the order of 20 % appear at M ~ 7.. 8 GeV (assuming the EKS98 
modifications) but below the J/\E' peak they are only about 5%. The data points at 
large masses, however, have a smaller weight in the x^ fits [J18l due to relatively large 
statistical uncertainties. The fits are dominated by pair masses near 4 GeV where the 
error bars are smaller. The nuclear effects in the mass distributions thus remain less 
than 5% for the extrapolation of the DY cross sections from 4 GeV down to 3 GeV. On 
the other hand, the experimental i^-factor, K = a^p/cr^j^ ' |T^ includes the region of 
large M, and thus could be underestimated by 10. .15 % if the nuclear effects in parton 
distributions are neglected. 
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